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Abstract — Reminiscent of the parity function in network coding 
for the butterfly network, it is shown that forwarding an even/odd 
indicator bit for a scalar quantization of a relay observation 
recovers 1 bit of information at the two destinations in a 
noiseless interference channel where interference is treated as 
noise. Based on this observation, a coding strategy is proposed 
to improve the rate of both users at the same time using a relay 
node in a noisy interference channel. In this strategy, the relay 
observes a linear combination of signals sent by the two sources, 
and broadcasts a common message to the two destinations 
over a shared out-of-band link of constant rate Ro bits per 
channel use. The relay message consists of the bin index of a 
structured binning scheme obtained from a 2 fl °-way partition 
of the squared lattice in the complex plane. We show that 
such scalar quantization-binning relay strategy asymptotically 
achieves the cut-set bound in an interference channel with a 
common out-of-band relay link of limited rate, improving the 
sum rate by two bits for every bit relayed, asymptotically at 
high signal to noise ratios (SNR) and when interference is 
treated as noise. We then use low-density parity-check (LDPC) 
codes along with bit-interleaved coded-modulation (BICM) as a 
practical coding scheme for the proposed strategy. We consider 
matched and mismatched scenarios, depending on whether the 
input alphabet of the interference signal is known or unknown 
to the decoder, respectively. For the matched scenario, we show 
the proposed strategy results in significant gains in SNR. For 
the mismatched scenario, we show that the proposed strategy 
results in rate improvements that, without the relay, cannot 
be achieved by merely increasing transmit powers. Finally, we 
use generalized mutual information analysis to characterize the 
theoretical performance of the mismatched scenario and validate 
our simulation results. 



I. Introduction 

Wireless communication is changing gears. The most 
promising way to meet higher and higher demands is to 
move towards a much higher spatial reuse, with very small 
cells where the distance between transmitters and receivers is 
dramatically reduced. On the other hand, this approach calls 
for a more and more unregulated user-based deployment of 
wireless networks, in contrast to the carefully and centrally 
planned conventional macro-cellular layout. In this context, 
femtocells (i.e., small home-based base stations) is rapidly 
emerging. Soon, people could be operating their own base 
stations (BS) in home or office, selectively controlling user 
access to their privately-operated BS. 
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Fig. 1 . Aside from the traditional role of data retransmission, relays can be 
used as "helpers" to resolve interference. In this figure, neighboring mobile 
stations (MS) in residential cellular communication networks inevitably inter- 
fere in the uplink, particularly because of restricted access policy in residential 
cellular systems (users cannot be distributed efficiently among BSs) and since, 
the stationary base stations cannot move to find a less-interfered spot. A relay 
station (RS) with a digital link to both base stations (BS) via Internet is used 
to reduce interference. See also a similar scenario in Fig. ^4] where the relay 
and base stations are integrated. 



Yet, interference remains as a major obstacle to widely de- 
ploy ad-hoc user-managed cellular systems, requiring further 
attention and new solutions. In a dense area, like a high- 
rise residential building, inter-cell interference is inevitable 
among user-controlled BSs. Because of sizable number of 
neighboring cells and small number of users in each cell, tradi- 
tional frequency reuse strategies are ineffective and spectrally 
inefficient. Hence, new strategies to handle interference are 
needed. 

This paper introduces a simple such strategy based on 
relays. Relays are traditionally used as middle terminals to 
retransmit the message. In an interference scenario, however, 
we envision an alternative application of relays as "helpers" to 
resolve interference. Fig. [T] shows an out-of-band relay placed 
between two neighboring interfering cells. The relay has a 
noiseless link of a limited rate to the two destinations, and can 
broadcast a common signal to the two destinations. The relay 
could be connected over an orthogonal wireless connection 
to the two users in the case of downlink (e.g., WiFi) , or 
even wired to the two BSs for the uplink communication (e.g., 
Ethernet). 

We focus on relay strategies that involve no decoding of 
the user message. This is because the relay also experiences 
interference, rendering the decoding difficult. In our strategy, 
the relay operation involves a symbol-level quantization of 
the relay observation using a planar lattice. The relay sends 
a bin index for its quantized observation, obtained through a 
suitable partition of the square lattice (T). See Fig.|4]and Fig. [5] 
for examples. The relay message is incorporated in the user 
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decoder to enhance the channel (initial) log-likelihood ratios 
(LLR) values. In other words, only initial LLR computations 
in the receiver decoder has to be modified to integrate the relay 
message in decoder. 

We simulate the performance of the proposed strategy using 
bit-interleaved coded modulation (BICM) and low-density 
parity-check codes (LDPC), assuming perfect channel state 
information at the decoder. In practice, channel parameters 
can be estimated using pilot subchannels with relatively small 
overhead, specially in femtocell scenarios with slow channel 
variations. We study matched and mismatched decoding met- 
rics, where the decoder is aware of the interference signal 
alphabet (constellation), or treats interference as Gaussian, 
respectively. In the matched scenario, where the decoder 
searches over the product of user constellations to compute 
LLRs, our results indicate that an SNR gain of approximately 
1 dB can be achieved for every bit relayed. The mismatched 
scenario is more interesting from a practical perspective, since 
it leads to a significantly lower decoding complexity, and our 
results indicate that a signficant gain is obtained using the 



relay. See Fig. 12 



The proposed relay strategy can be viewed as an analog 
version of digital network coding for wireless channels |2j, |3J. 
We explain this connection to network coding and the relation 
of the proposed scheme with respect to the compress-and- 
forward (CF) relaying scheme (4|-|(6| in Section [II] We show 
that for an interference channel where interference is treated as 
noise, a scalar quantization strategy is asymptotically optimal. 
That is, we show that asymptotically in the low noise regime, 
every bit relayed improves the achievable sum rate by two bits, 
using a scalar quantize-and-forward strategy. This two-for-one 
improvement in sum rate is analogous to the rate improvement 
due to parity forwarding in the celebrated butterfly network 
example in network coding |7). 

The interference channel with an out-of-band relay has 
been the subject of a number of recent studies, where the 
fundamental information theoretic aspects of the interference 
relay channel are investigated p)-fTT). However, practical 
coding strategies to utilize relays in interference scenarios 
are not as widely studied. Moreover, not many practical 
compress-and-forward relay coding schemes are available, 
since the decode-and-forward scheme is often favored upon 
the compress-and-forward relay strategy for multihop relaying. 
In this work, we make further progress in bringing some of 
the theoretical achievements on the application of relays in 
interference scenarios to practice by devising a simple, yet 
effective, coding strategy based on the compress-and-forward 
scheme, tailored specifically for the interference channel. 

The rest of this paper is organized as follows: Section [TT] 
describes some of the connections between the proposed 
strategy, digital network coding, and compress-and-forward 
relay strategy. Section [III] illustrates the system model and 




Fig. 2. A two-user interference channel with a common out-of-band relay 
link of rate Ro. 
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Fig. 3. Coarse parity: the relay forwards a parity bit for a quantized value 
ofY r . 



encoding and decoding procedures, and Section VI provides 
some simulation results. Finally, a few remarks conclude the 



paper in Section VII 



II. Background 

A. Network Coding 

A simple illustration of coarse network coding can be given 
using a toy example. Consider a noiseless version of the 
interference channel shown in Fig. [2] with N\ — N 2 = 
N r = 0. In this model, the source signal Xi is corrupted 
with interference X 2 and is received as huXi + h 2 \X 2 at 
the destination, with X 2 treated as noise. Assume that all 
channel gains are available at destination nodes and Xi and 
X2 are real independent Gaussian random variables. Consider 
now a relay that observes Y r = g\X\ + g 2 X 2 and wishes to 
assist the destination by forwarding 1 bit of information over 
a shared noiseless link. A simple relay strategy to improve the 
achievable rate at the destination by close to one bit is to send 
a if L^r /rfJ is even, or a 1 otherwise, for a very small value 
of d > 0; see Fig. [3] To see this, let X r = [Y r /d\ mod 2 
denote the relay message, and then the rate improvement is 
given by: 

I(X 1 ;Y 1 ,X r )=I{X 1 ;Y 1 )+I(X 1 ;X r \Y 1 ) 

= J(Xi;Ki) + I(Xx;X r \huXi + h 21 X 2 ) 
= J(Xi;Ki) + H(X r \Y 1 ) - H{X r \X 1 ,Y l ) 

( =' /(Xi; Yi) + HiXrlhnX, + h 21 X 2 ) - 

( = } /(X i; Yi) + l 

In the above derivation, (a) follows since Y r is known given 
Xi and Y\ and thus, H(X r \X\, Y\) = 0, and (b) follows since 
for very small values of d, \Y r /d\ mod 2 is a Bernoulli 1/2 
random variable given Y\, provided that the relay observation 
is not statistically identical with the receiver observations, i.e., 
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the matrix 

hn gi 

is full rank. (See Appendix.) 

Now, switch the roles of Xi and X 2 , and consider the 
second destination who is interested in decoding X2, while X\ 
is now interference and is treated as noise. Interestingly, the 
same relay strategy improves the achievable rate of the second 
user also by 1 bit. In other words, a single bit from the relay 
recovers one bit of information at each of the two destinations. 
This is an example of coarse network coding, where an 
even/odd indicator for a quantized version of relay observation 
recovers two bits of information, and is reminiscent of the 
celebrated butterfly network example in network coding (7). 
A simple generalization of the above coarse parity function 
can be used to allow for higher relay data rates, and more 
importantly, to account for background noise. 

Although focusing on different problems, there are in- 
teresting connections between the proposed coarse network 
coding and the analog network coding strategy devised in fl2) . 
Analog network coding proposed in [ 12] is a symbol-wise 
nonlinear amplify-and-forward, where a scalar relay function 
is optimized (assuming differentiability) such that the end-to- 
end mutual information between the source and destination 
is maximized. Interestingly, the optimized nonlinear relay 
functions have a semi-periodic form, resembling a smooth 
version of the scalar quantization and binning of coarse net- 
work coding (see Section |III- A[ >. The analog network coding 
strategy, however, is devised only to serve a single destination 
since the relay operation depends on the destination channel, 
and the system optimization searches only for differentiable 
relay functions, which may not be necessarily optimal. 

The problem considered in this paper differs from the 
wireless network coding approach of (2), where it is assumed 
that the relay observes clean versions of the source data. 
Coarse network coding also differs from the compute-and- 
forward scheme of [13] in that structured codes at the sources 
are not required in coarse network coding as no decoding, even 
of a function of the two source messages, is performed at the 
relay. Further, coarse network coding is a symbol-wise strategy 
in the same spirit of the XOR strategy of digital network 
coding. 

B. Compress-and-Forward Relaying 

Coarse network coding is essentially a scalar compress- 
and-forward strategy ( p] Theorem 6]) where quantization 
and binning are performed at the symbol level. Though a 
scalar quantization and binning are not generally the strategies 
of choice for a single-relay channel, this paper shows that 
there are substantial gains to be obtained from a simple scalar 
compress-and-forward strategy in an interference channel. 

The asymptotic incremental optimality of scalar compress- 
and-forward shown in the previous section (i.e., one bit rate 
improvement per one bit relayed as noise tends to zero) is a 
consequence of cross-determinism in a noiseless interference 
channel. A cross deterministic relay channel was first intro- 
duced in p4| , and consists of a three-terminal network of a 



source X, a relay Y r , and a destination Y, with a deterministic 
relationship between the relay observation Y r , and the sent and 
received signals X and Y, in way that Y r is a deterministic 
function of X and Y. The noiseless interference channel in 
Section is an instance of two coupled cross-deterministic 
relay channels: here, from X\ and Y\ or from X2 and Y 2 , we 
can compute the relay observation Y r ; for both users, the relay 
channel is cross-deterministic. 

To see the role of cross-determinism in the compress-and- 
forward strategy consider the achievable rate of compress-and- 
forward with list decoding for the interference channel with 
interference treated as noise |6): 

Ri < I{Xi ; Y x ) + min{i? , 1(Y r ; Y^)} - I(Y r ; Y r \X U Y t ) 
R2 < I{X 2 -Y 2 ) + mia{R ,I{Y r ; Y r \Y 2 )} - I(Y r ; Y r \X 2 , Y 2 ), 

(1) 

where Y r is an auxiliary random variable that plays the role 
of quantization of the relay observation Y r . When Y r is a 
deterministic function of X\ and Y\, and also X 2 and Y 2 , the 
penalty terms in ([T} vanish, since 

I(Y r ;Y r \X 1 ,Y 1 ) = 
I(Y r ;Y r \X 2 ,Y 2 ) = 0. 

Consequently, the rate of each user is improved by Rq bits, as 
long as the minimums in ([l) occur at Rq for the quantization 
scheme used. In other words, performance is not very sensitive 
with respect to the quantization scheme and in particular, the 
block length over which quantization is performed. 

Unfortunately, optimal vector quantization with side in- 
formation and list decoding as developed in [5] or [6| is 
not quite amenable to practical code construction. The list 
decoding scheme of (6}, or the more recent noisy network 
coding approach of [5], search over all possible quantized relay 
observations Y r and transmitted source codewords over a large 
block. In other words, an exponential number of codewords 
have to be examined in order to decode the source message. 

However, this paper shows that much of the promised 
gain can be achieved at high SNRs by simply using scalar 
quantization and binning. This is because the requirements 
on Y r to achieve i?o bits of improvement are quite loose at 
high SNRs, in the sense that Y r is admissible as long as it is 
sufficiently close to Y r (measured in terms of I(Y r ;Y r \Yi) 
and I(Y r ;Y r \Y 2 )). Driven by the above intuition, we can 
replace the high-dimensional vector quantizer of the general 
CF scheme with a simple scalar quantizer. Most remarkably, 
when this strategy is combined with conventional (single- 
user) powerful coded-modulation based on binary LDPC with 
BICM and iterative decoding, and the iterative decoder inputs 
are enhanced by incorporating the relayed bits, it is shown 
that a significant gain is achieved even at practical (non- 
asymptotically high) SNR. 
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III. System Design 

Fig. [2] shows the baseband-equivalent model of a two-user 
interference channel with a relay defined as: 



Y\ — huXx + h 2 \X 2 + N\ 

Y 2 = h 22 X 2 + ^ 12-^1 + N 2 , 

with a relay observing 

Y r = gi Xi + g 2 X 2 + N r , 



(2a) 
(2b) 



(2c) 



where Xi and X 2 are source signals of maximum power Pi 
and P 2 , and Y\, Y 2 , and Y r denote the channel observation 
at the two destinations and the relay, respectively. The relay 
assists both destinations simultaneously by broadcasting X r 
over a common digital link of rate Rq. Here, N%, N 2 , and N r 
are circularly-symmetric complex Gaussian noise of zero mean 
and variance Nq/ '2 per real dimension, and hij,gi,i,j = 1,2 
represent channel gains. We assume that the channel undergoes 
fast fading with independent channel gains across different 
time slots, and hij,gi,i,j — 1,2 are modeled as independent 
circularly-symmetric complex Gaussian random variables of 
zero mean and variance 0.5 per real dimension. Also, rele- 
vant channel state information (CSI) is known to each user, 
specifically, g\ and g 2 are known to the destinations and the 
relay , and h 22 , h\ 2 are known to destination 1 , and h 22 , 
h\ 2 are known to destination 2 . In practice, the i.i.d. fast 
fading model with CSI known to all can be approached by 
using interleaving and pilot-assisted channel estimation, for 
a time- and frequency-selective fading channel with OFDM 1 , 
where sufficiently long codewords are interleaved in the time- 
frequency domain, as currently done in today's OFDM sys- 
tems. 
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Fig. 4. Chessboard quantization of Y T at rate Rq = 1. 
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Fig. 5. Chessboard quantization of Y r at rate Rq = 2. 



A. Encoding: 

The source message is encoded using a conventional LDPC 
code, then the binary codeword is interleaved and mapped to a 
sequence of M-ary constellation symbols using Gray labeling. 
The sequence of constellation symbols is sent over the channel 
using quadrature amplitude modulation (QAM). In ©, X\ and 
X 2 represent the constellation symbol sent by the two users 
in each time slot, and X r represent the relay message of rate 
R , i.e., H{X r ) = Rq. 

In each slot, X r is computed according to a scalar quantiza- 
tion scheme that we refer to as chessboard quantization, since 
the relay quantization Voronoi pattern resembles a chessboard. 
Chessboard quantization can be described as coset partitioning 
of the square lattice into 2 Ra sublattices [1|. Denote the the 
square integer lattice in the complex plane by Z 2 , and let A 
denote a nested sublattice of Z 2 with generator matrix G and 
nesting ratio of 2 R °, i.e., 

A = {G-i:ieZ 2 }, (3) 

where G is a matrix of integers with 

det(G) = 2 fl °. (4) 

'Orthogonal Frequency-Division Multiplexing. 



We use A to partition Z 2 . Chessboard quantization can then 
be formally defined as: 

X r = B [Q(Y r /d)\ , (5) 

where d is a complex scaling factor, Q(-) represents the nearest 
neighbor quantizer with respect to Z 2 , and B[ i ] denotes the 
coset leader for i € Z 2 with respect to the coarse lattice A. 

For given i?o, gi,g 2 , Nq, and input alphabets, we need to 
optimize the quantization, i.e, d, and A. For practical feasibil- 
ity, we restricted the quantization lattice to be Z 2 and we have 
optimized A independently of the channel realization, i.e., A is 
chosen depending on the statistics of the channel coefficients 
and kept fixed for all channel uses. A more complex alternative 
would be to choose A among all sublattices with nesting ratio 
2 R " based on the value of the channel coefficients gi,g 2 . 
However, this more complex solution was not pursued here. 
In summary, only the scaling coefficient d is optimized on a 
symbol-by-symbol basis, based on the realization of gi,g 2 . 
In principle, the best partition for given source constellations 
and channel gains may be very complicated if we remove 
the lattice structure. Fixing the sublattice partition pattern first 
significantly reduces the optimization complexity. 

To find a suitable lattice partition, we can optimize A such 
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that the normalized second moment of A is maximized [ 1 1 . For 
special cases Rq = 1, 2, we choose A as the lattice generated 
by 



for Rq = 1, and 



G = 



G 



(6) 



(7) 



for R = 2. See Fig. g] and Fig. [5] 

The parameter d is found such that H(X r \Xi, X 2 ) for given 
gi,92, No and input alphabets is minimized. An explanation 
for choosing H(X r \Xi, X 2 ) as the optimization metric is 
given later based on log-likelihood ratios. As an alternative 
motivation, minimizing H(X r \Xi, X 2 ) is related to maximiz- 
ing the achievable rates. For a given channel realization, the 
rate improvement due to the relay for user one and user two 
is given by: 

Ai?i = I(Xi; X r \Yi) = H{X r \Yx) - H{X r \X u Y x ), 
AR 2 = I(X 2 ;X r \Y 2 ) = H(X r \Y 2 ) - H(X r \X 2 ,Y 2 ). 

Since the relay sends a common signal to both destinations 
and direct channel gains hij,i,j = 1,2 are not known at 
the relay, Ai?i and AR 2 cannot be jointly maximized at the 
relay for each time slot. An alternative would be to minimize 
H(X r \Xi, X 2 ) for given g\ and g 2 (which are available at the 
relay), since for a given channel realization: 



HiXrlXuYi) > H(X r \X 1 ,X 2 ,Y 1 ] 
H(X r \X 2 ,Y 2 ) > H(X r \X 1 ,X 2 ,Y 2 ) 



(b) 



: H(X r \X 1 ,X 2 ) 
H(X r \X u X 2 ) 



where (a) follows from conditional entropy inequality, and (b) 
follows since for given channel gains, X r — (Xi,X 2 ) — Y\ 
forms a Markov chain. Thus, H{X r \X\, X 2 ) is a lower bound 
on both H(X r \X 1 ,Y 1 ) and H(X r \X 2 ,Y 2 ), which becomes 
tight at high SNR where background noise is negligible. Con- 
sequently, minimizing H(X r \X\,X 2 ) asymptotically results 
in maximized Ai?i and Ai?2- 

IV. Decoding 

This section describes the decoding scheme for user one. 
Decoding at user two follows similar steps. Given Yj and 
X r , destination one computes LLRs for corresponding bit 
positions of the underlying LDPC code, and then conventional 
sum-product algorithm is used to iteratively decode the source 
binary codeword. In other words, X r only affects the initial 
computation of LLRs at the destination. 

We consider two forms of decoding: A matched scheme 
where the actual statistics for the interference signal is used 
for decoding, and a mismatched scheme where the distribution 
of the interference signal is approximated to be Gaussian 
with the same mean and variance. For QAM modulation, the 
matched decoder searches over the Cartesian product of the 
source and interference constellations to compute LLRs. Such 
an exhaustive search over the product constellation generally 
achieves a better performance at the cost of higher complexity. 



On the other hand, the mismatched decoder searches only over 
the source constellation and approximates the LLRs by treating 
the interference signal (of discrete-alphabet) to be Gaussian. 

In practical settings, the interference signal constellation is 
often unknown to the decoder, since the control channel of 
the interfering user usually needs to be decoded to obtain the 
encoding parameters like the constellation size. In such sce- 
narios, the mismatch decoding scheme is an attractive strategy 
with a lower decoding complexity and robust performance. 

A. Matched Decoder 

Let bi ■ ■ ■ 6fc denote the binary Gray label for constellation 
symbol s € A4, where M. is a constellation of size M = 2 . 
Define also the reverse mapping 6j = Ai(s),i — l,...,k, 
where 6j denotes the z'th bit of the binary label of s. 

Given yi and x r , the destination computes A*, the LLR 
corresponding to the z'th bit position, for i = l,...,k as 
follows: 

^ _ j Y,{x 1 £M 1 ,x 2 £M 2 \A l (x 1 ) = l} P(Vli X r\ x l, %2) ^ 
^2{x' 1 EM u x' 2 <EM 2 \A,{x' 1 )=0} P(V^> X r\%l, X 2 ) 

where A4i,M 2 are the constellations for user one and two. 
Now, to compute p(yi,x r \xi,x 2 ), we have: 

p(yi,x r \xi,x 2 ) = p(yi\xi,x 2 )p(x r \xi,x 2 ), (9) 

where 

1 f 1 1 2/1 - hnxi - h 2 ix 2 \\ 2 ' 



p(yi\xi,x 2 ) 



itNq 



exp 



and 



1x1,22) 



ttN, 





y r £B(x r ) 



cxp 



N 



91x1 - g 2 x 2 \ 



dy r 



where B(x r ) encompasses all quantization regions with X r = 
x r ; for example, B(0) represents all gray-colored squares 
corresponding to x r = in Fig. [4] 

The LLR functions in |8]l and (|9]l also suggest that minimiz- 
ing H(X r \Xi, X 2 ) amounts to better initial LLRs with larger 
magnitude. This can be illustrated by where the effect of 
the relay signal X r on the destination LLRs is controlled by 
p{x r \x\, x 2 ). Hence, the more deterministic p(x r \xi,x 2 ) is, 
the larger the magnitude of initial LLRs are, which results in 
a lower decoding error rate. 

B. Mismatched Decoder 

Computing the LLR values in ([8]) requires summing over 
all pairs of symbols in and M 2 , which can be com- 
putationally intensive depending on the size of .Mi and M. 2 . 
This section describes a low-complexity mismatched decoding 
scheme that does not require the knowledge of M. 2 . 

The mismatched decoder treats the interference signal X 2 as 
a Gaussian random variable with the same mean and variance. 
Yet, the mismatched scheme becomes matched when X 2 is 
indeed Gaussian; for example, asynchronous cross sub-channel 
interference in an OFDM system due to misaligned timing can 
be well approximated to be Gaussian. 



6 



Decoding is again performed by computing the initial 
mismatched log-likelihood ratios (mLLR). Given y\ and x r , 
user one computes 7$, the mLLR corresponding to the i'th bit 
position, for i = 1, . . . , k as follows: 

Y,{ Xl eM 1 \A i (x 1 )=l} q{yi, x r\xi) 

7i = lo g / fTT' ( 10 ) 

2j{x' 1 eA'fi|Ai(x' 1 )=0} 1\V\i x r\ x i) 

where Ai 1 is again the constellations for user one. The metric 

^(yijavlxi) is given by: 



q(yi,x r \xi) = q\(y\\x{)qr(x r \xi, yi), 



(11) 



where 

qi(yi\x\) 

and 



^(||/l2l|| 2 P2+^0) 



exp 



\\yx - hnxi\\ 2 
\h 21 \\ 2 P 2 +N 



g r (a; r |a;i,yi) 



exp 



\Vr - Mil 5 



rfy r , (12) 



y,,£B(x 



with 



a z =N + 

/! = 51X1 



\\92fP2No 

\h 21 \\ 2 P 2 + N 

92h* 2l P2 

' \\h 2 i\\ 2 P 2 + N, 



(yi - huxij. 



Note that q(yi, x r \x r ) — p(yi, x r \x\) when X 2 is Gaussian. 
The mLLRs computed for every bit position are used as initial 
LLRs fed to the receiver's iterative decoder. 

V. Achievable Rates 

Under the matched decoding scenario, the achievable rates 
of the above strategy can be computed using standard mutual 
information analysis. However, characterizing the capacity 
limits under the mismatched decoding scenario requires a 
different set of techniques. Based on the analysis in this 
section, numeric results for the achievable rates under matched 
and mismatched scenarios are presented in Section |VT] 

A. Matched Metric 

In the matched scenario, the decoder uses p{y il x r \xi) as 
the likelihood metric, i = 1, 2 for user one and two. For this 
decoding metric, the achievable rate is given by the mutual 
information between channel input Xi and the channel outputs 
Yi, X r for i = 1,2. Thus, we have 

p(X 1 ,Y 1 ,X r ) 



i?i =I(X 1 -Y 1 X r ) = E 
R 2 =I(X 2 ;Y 2 X r ) = E 



loe 



log 



piXMYuXr) 
p(X 2 ,Y 2 ,X r ) 
p(X 2 )p(Y 2 ,X r ) 



(13) 
(14) 



where the expectation E is also with respect to channel 
coefficients 2 . Using ( |13) , R\ can be computed for discrete- 

2 A11 the probability distributions and mutual information 
expressions in this section are conditioned with respect to the 
channel coefficients, but explicit conditioning is dropped for 
notational simplicity. For example, p(yi,x r \xi) should be treated as 

PY i ,X r \X i ,h 1 i,~;9i,92(Vi' Xr \ x i> hl1 ' ■ ■ -^22,91,92) and I(X 1 ;Y 1 X r ) 



alphabet X\ and X 2 inputs using 

p(y 1 ,x 1 ,x r )= ^2 P(3Jii x r\xi,x 2 )p{xi)p(x 2 ), (15) 

along with ([9j, where p(xi) and p(x 2 ) are uniform over the 
source constellations. When X 2 is continuous, an equivalent 
expression can be found for p(yi, x\, x r ) by using integral in 
place of sum in ( [15) . The rate R 2 can be computed similarly. 

B. Mismatched Metric 

The mismatched decoder uses q(yi, x r \xi) as the likelihood 
metric. The supremum of achievable rates of a mismatched 
decoder is unknown in general. However, error exponent 
analysis with random coding can be used to establish an 
achievable rate for a mismatched decoder p5) . 

Theorem 1 (Generalized Mutual Information ^15^ ): 
Consider a point-to-point channel defined by p(y\x) for 
X 6 X and Y 6 y. Using a mismatched decoder with 
mismatched likelihood metric q(x,y), all positive rates R are 
achievable for this channel with a vanishing error probability 
for large block lengths if 



R < I am, (X; Y) ^ max I 9 s mi {X; Y). (16) 

s>0 



where 



jgrm A £ 



log 



q(X,YY 



E x <exP( x ')l( x ',Yy 



(17) 



Using Theorem [T] computing lower bounds for the achievable 
capacity with the mismatched decoding metric in ( fTT| is quite 
straightforward. For user one, we have 



i?i(s) = E 



lop 



q(Yi,x r \Xi) 6 



(18) 



denotes E 



h 11 ,...,h 2 2,91,92 



I(X 1 ;Y 1 X r \h 11 ,...,h 22 ,g 1 ,g 2 ). 



where Mi is the size of A4i. The expectation E is with respect 
to probability distribution p(yi, xi,x r ) given in ( |15) , and also 
the channel coefficients. To compute i?i(s) numerically with 
fading coefficients, the average I§ mi is computed for different 
realizations of the channel coefficients, with s optimized for 
each realization. A similar expression can be found for user 
two. 

VI. Simulation Results 

This section first investigates the theoretical achievable 
rate for the matched and mismatched scenarios. Then, the 
performance of iterative decoding for LDPC codes with BICM 
modulation is studied using simulations. 

A. Achievable Rates 

For the matched decoding scenario, Table [I] lists the im- 
provement in minimum SNR required to achieve a source rate 
of 1 bit/symbol for different constellation sizes as channel 
inputs. As it is shown, when both the desired signal and in- 
terference are taken from 4-QAM constellations, an SNR gain 
of 0.86 dB is obtained with 1 bit relayed. With 2 bits relayed, 
the SNR gain is more substantial: 2.75 dB improvement to 
achieve a data rate of 1 bit/symbol. The gain due to coarse 
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TABLE I 

Improvement in Minimum required SNR to 
achieve ill = 1.5. 



Mi 


M 2 


Ro = 1 


R = 2 


4-QAM 


4-QAM 


0.86 dB 


2.75dB 


4-QAM 


16-QAM 


2.03 dB 


4.74 dB 


16-QAM 


16-QAM 


2.18 dB 


4.78 dB 




■ No relay 

■ 1 bit relayd 
2 bit relayed 



30 



15 20 
SNR(db) 

Fig. 6. Achievable rate with matched decoding and QPSK Xi and X2. 



network coding improves at higher bit-per-symbol rates. The 
second and third row of Table [IJlist the gain in minimum SNR 
required to achieve R\ = 1, when interference X2 is taken 
from a 16-QAM constellation. In this case, 1 bit of relaying 
results in about 2 dB gain in SNR, and the SNR gain is almost 
doubled in dB with 2 bits relayed. 

Fig. [6] shows the achievable rate i?i when both X\ and X2 
are taken from a 4-QAM (QPSK) constellation of power P, 



_ 1 6QAM/16QAM, R =0 
_ 1 6QAM/16QAM, FL.1 
1 6QAM/16QAM, R .2 

- 4QAM/16QAM, R o =0 

- 4QAM/16QAM,R .1 
4QAM/16QAM.R..2 




T 16QAM (R Q =0; 

64QAM fR Q =0; 

, 16QAM (R Q =1 

— j 64QAM (R Q =1 

16QAM (R Q =2; 
□ 64QAM (R =2; 




SNR(db) 



Fig. 7. Achievable rate with matched decoding and 16-QAM X2. 



Fig. 8. Achievable rate with matched decoding for discrete-alphabet X\ 
and Gaussian interference X'2- 



and SNR is defined as 10\og 10 (P/N ). With discrete input 
alphabets, the destination can uniquely identify the source 
symbol sent by X\ at high SNRs, and thus asymptotically, the 
achievable rate Ri tends to log 2 (4) = 2, as shown in Fig. [6] 
Fig. [6] also indicates that a substantial gain in minimum SNR 
required to achieve a certain rate is achieved. Highest SNR 
gains are obtained for larger values of R\. At SNR=25 dB, 
the SNR gain is close to 4 dB for one bit relayed, and 8 dB 
for 2 bits relayed. 

Fig. [7] shows the achievable rate curves when interference 
comes from a 16-QAM constellation. In this figure, X\ and 
X2 are of the same power P, and SNR is 10 \og 10 (P/No). 
Comparing the first and second row in Table [IJ it is revealed 
that coarse network coding is more helpful for a QPSK source 
when interference X2 comes from a 16-QAM constellation as 
opposed to 4-QAM. A possible conclusion is that the more 
damaging the interfering signal is, the more the relay can help. 

Fig. [8] shows the rate improvement obtained by coarse 
network coding when the source signals X\ is taken from 
a 16-QAM and 64-QAM constellations and interference is 
a circularly-symmetric complex Gaussian random variable. 
Both Xi and X 2 are of the same power P and SNR is 
again 10 \og 10 (P/N ). In this case, coarse network coding 
almost linearly improves the achievable rate. Highest gains 
are achieved when X\ is taken from a 64-QAM constellation, 
where asymptotically at high SNRs every relayed bit improves 
the achievable rate by slightly less than 1 bit. 

From Fig. [6] Fig. [7] and Table [I] it is expected that coarse 
network coding improves the achievable rates by larger extents 
as the size of the source alphabets increase. Fig. [9] shows 
the achievable rates when both X\ and X2 are circularly- 
symmetric Gaussian random variables. It is shown in Fig. [9] 
that for Gaussian inputs, the achievable rate of both users is 
improved by close to 1 bits at Rq = 1, and 1.5 bits at Rq = 2, 
at high SNRs. 



Fig. 10 shows the achievable rate obtained by the general- 
ized mutual information analysis for the mismatched decoding 
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Fig. 9. Achievable rate with matched decoding for Gaussian inputs. 
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Fig. 10. Achievable rate with mismatched decoding and QPSK inputs. 



metric and QPSK inputs. Although both the inputs X\ and X 2 
are discrete, the achievable rate with mismatched decoding is 
saturated even at high SNR values due to interference. Yet, 
adding the relay alleviates the rate saturation at high SNRs by 
improving the asymptotic achievable rate in proportion to the 
number of bits relayed. The improvement in achievable rates 
for QPSK inputs and mismatched decoding is less than a bit for 
every bit relayed. However, this rate improvement tends to one 
as the size of signal constellations increases. At the extreme 
case where both X\ and X 2 are continuous Gaussian random 
variables, the mismatched decoding reduces to the matched 
strategy, and every bit relayed results in two bits improvement 
in sum rate for large SNR. 

B. Performance with LDPC-Encoded BICM 



For the matched scenario, Fig. 11 shows the performance 
of BICM-LDPC coding along with coarse network coding 
and iterative decoding at source rate R% = 1. The BICM 



SNR(dB) 

Fig. 11. The bit-error-rate performance of regular (3,6) LDPC code with 
matched decoding and BICM gray-mapped signaling. Block length is 20000. 



system is comprised of an LDPC code of rate 0.5 and block 
length 20,000, and gray-labeled 4-QAM modulation for both 
X\, and interference X%. The underlying LDPC code is 
the Gallager regular (3,6) LDPC code flT6}, with a random 
graph construction. Decoding is performed as described in 
Section [IV] by enhancing the initial LLRs using the extra bit 
received from the relay. 

An SNR gain of approximately 1 dB is achieved using the 
(3,6) LDPC code of rate 1 /2 along with coarse network coding 
at R = 1. At R = 2, the gain in SNR is approximately close 
to 2.5 dB. Improved performance could perhaps be achieved 
by optimizing LDPC codes for the specific enhanced LLR 
density distribution. Because of the decoding approach used, 
the achieved gain is controlled by the shape of enhanced 
LLR density distribution. The more the relay message skews 
the initial LLR density distribution, the higher the achieved 
SNR gain would be, which also requires tuning the LDPC 
code for the specific input LLR density distribution to fully 
exploit potential gains. In this paper, standard LDPC codes 
are used to simulate the performance, as our focus is more 
on the multiuser aspects, rather than the impact of LDPC 
graph itself. Also, since the relay message is only used for 
initial LLR computations, the iterative decoding procedure at 
each destination is equivalent to that of a single-user point- 
to-point channel, which can be optimized by standard tools 
such as density evolution, EXIT charts, and protograph-based 
structures p7l-fT9l. 

Fig. [12] shows the performance of the (3,6) LDPC code of 
rate 0.5 with mismatched decoding and gray-labeled QPSK 
BICM signaling. Using the relay results in substantial im- 
provements when the decoder is unaware of the interfering 
signal constellation. At rate Ri = 1, receiver is unable to 
decode the source codeword with mismatched LLR values, 
even at high SNRs, without relay assistance. Yet, using a 
relay link of rate Rq = 1 enhances the LLR values in a 
way that decoding is possible at SNR values around 10 dB, 
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SNR(dB) 

Fig. 12. The bit-error-rate performance of regular (3,6) LDPC code with 
mismatched decoding and BICM gray-mapped signaling. Block length is 
20000. 



which is consistent with GMI achievable rates in Fig. 11. 
A relay link of rate Rq = 2 results in an additional 5 dB 
SNR gain as compared to i?o = 1. again matching the result 
from GMI simulation in Fig. 



10 The SNR gains are quite 



substantial, especially if we consider the very simple relay 
strategy, requiring no decoding and only scalar quantization. 

C. An intuitive interpretation 

It is helpful to intuitively understand how ([8]) and ( fTTj ) 
enhance the LLR values at the decoder. We can gain some 
insights by looking at likelihood computation for a particular 
constellation symbol Sj sent by user one. Assuming that 
si is sent by user one, the decoder can subtract si from 
its observation y% to find an estimate for S2, the symbol 
sent by user two. Now, using the channel state information, 
decoder can find an estimate for the relay observation and 
the corresponding quantization index. If the decoders estimate 
of the relay message matches the message received from the 
relay, the likelihood of s% being sent is improved, otherwise, it 
is less likely that si is sent by user one. Such small adjustments 
to the LLR values at symbol level can be combined over a 
larger block length using, for example, an LDPC code with 
iterative decoding. Small enhancements in initial LLR values 
translate to significant improvement in overall bit error rate 
(BER) performance, because of the waterfall characteristic of 
BER curves for LDPC codes. 

VII. Conclusion 

Interference is one of the main obstacles to deploy user- 
based uncoordinated wireless cellular networks. We propose 
a simple relay strategy to improve the achievable rates in a 
two-user interference channel, where interference is treated 
as noise. This relay strategy is called coarse network coding, 
since the relay operation involves a scalar quantization fol- 
lowed by a scalar binning, reminiscent of the parity generation 
in digital network coding. In addition, asymptotically at high 




Fig. 13. Each user can play the role of a relay for another user. The same 
chessboard quantization strategy can be used for limited receiver cooperation 
to resolve interference. 



SNRs, a single bit sent by the relay in coarse network coding 
improves the rate of both users by one bit, resembling the 
rate improvement obtained from a single parity bit at two 
different destinations in digital network coding. In this relay 
scheme, the relay quantizes its observation using a scalar 
quantization obtained by suitable partitioning of the square 
lattice. We consider decoding scenarios using matched and 
mismatched metrics, depending on whether where the alphabet 
of interfering signal is known or unknown to the decoder, 
respectively. Using LDPC-encoded BICM modulation, we 
show that the expected theoretical SNR gains due to coarse 
network coding are mostly achieved using a generic regular 
(3,6) LDPC code; higher SNR gains are expected by designing 
the LDPC graph for the specific density distribution of the 
relay-enhanced LLRs. This is left for future work. 

The proposed relay strategy can also be used to allow 
receiver cooperation in an interference channel. This is il- 
lustrated in Fig. [T3] where each receiver is connected via a 
digital link of limited rate to the other. Each receiver plays 
the role of a relay for the other user, quantizes its observation 
prior to decoding using the same chessboard-like quantization 
strategy, and forwards the index corresponding to its quantized 
observation to the other user. Decoding is performed again by 
computing enhanced LLR values and iterative decoding. 

Note that the analysis for the shared relay scenario shows 
that the same relay message can be used by different receivers; 
thus, generating relay messages are independent of the end 
receiver and the relay messages could be consolidated in 
a server to be used by any relevant interfered receiver. As 



an application, Fig. 14 shows a possible uplink scenario in 
a residential cellular wireless communication network where 
the two base stations use their Internet connections to com- 
municate with each other as limited-rate relay nodes via a 
server. In the uplink, the base station playing the role of 
a relay shares its channel state information and bin indices 
for its quantized observation with neighboring base stations 
using, for example, an interference resolution server. A base 
station experiencing interference queries the server for relayed 
information from neighboring interfering base stations. The 
receiver then incorporates the retrieved information to enhance 
the initial LLR calculations, which are then used by the 
iterative decoder. 
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Fig. 14. An uplink scenario where the two users cooperate as relays for 
each other to reduce interference via their Internet connection. Since the relay 
operation in coarse network coding strategy is independent of the destination 
being served, and since the same relay message can be used by different 
destinations, an interference resolution server (IRS) can consolidate relay 
messages from participating base station relays. An interfered receiver queries 
the IRS server to receive relevant relayed information about neighboring 
interference sources to assist its decoder. 



Appendix 

In this appendix, we prove the asymptotic incremental 
optimality of scalar quantization/binning for a Gaussian in- 
terference channel for Rq = 1, and with Gaussian inputs, 
where interference is treated as noise. That is, we prove that 
one bit relayed improves the rate of each user by one bit 
asymptotically at high SNRs. 

The incremental rate improvement for user one can be found 
as follows: 

I(X x ;Y x ,X r ) = I{X X ;Y X ) + I(X x ;X r \Y x ) 

= I{X X ; Yi) + I(X x ;X r \h xx X x + h 21 X 2 ) 
= I{X X -Y X ) + H(X r \Y x ) - H(X r \X 1 ,Y 1 ) 
= I(X X ;Y X ) + H(X r \h xx X x + h 21 X 2 + Ni) 
- H(X r \h lx X x + h 2X X 2 + N X ,X X ). 

We prove that as Nq — > 0, 

H(X r \h xl X x + h 2X X 2 + N x ) -)■ 1 (19) 



and 



H{X r \h xx X x + h 21 X 2 + N X ,X X ) -> 0. 



(20) 



Since X r = [Y r /d\ mod 2, we first need to find the 
conditional distributions of Y r given Y x and (Y X ,X X ). We 
have: 



r i s 1 J Ik-- Mill 2 ', , 

P{Vr \X X ,Vl) = exp <^ 2 f d Vr ' 

7T(T I erf 



(21) 



where 



and 



ai=N Q + 
Mi = 9%xi 

p(y r \yi) = 



\\92\\ 2 P2N 

\h 21 \\ 2 P 2 + N 

92h 2X P 2 
' \\h 21 \\ 2 P 2 +N 



(yi - h xl x x ) 



na 2 



exp 



\Vt - j"2 || S 



dy r 



(22) 



where 

2 M ||ffi^i -g 2 feiiH 2 PiP 2 +iVo(||.9i|| 2 Pi 
a 2 = 1\ Q H 



.1.92 



\ 2 P 2 ) 



V-2 = 

Let 



\\h xx \\ 2 P x + \\h 2X \\^P 2 + N Q 
g x h\ x Pi +g 2 h* 21 P 2 



\h lx \\ 2 P x 



\h 21 \\2P 2 



N 



yi- 



Internet 



with < a < 0.5. Then, we have: 



var(F r /d|Xi,Yi) 



0, 



!l 52 



\ 2 P2 



\h 2 i\\ 2 P2 



(23) 



as A^o — > for a < 1/2. Thus, given X x and Y x , Y r /d is 
asymptotically deterministic, and consequently, [Y/dJ and 
also [Yr- / dj mod 2 are deterministic. This proves (20) . 
To prove ( [T9) , notice that: 

12 i 



v&r(Y r /d IX^YA = 1 



\ 2 Pi + \\ 9 2rP2 



llffi 

|/lll|| 2 Pl + ||^2l|| 2 P2+A 

||gife 2 i-g 2 feii|| 2 -PiP2 

||/lll|| 2 Pl + ||/*2l|| 2 P2+A r 



• N, 



(i 



— > oo, 



■N 
(24) 



-2a 



as Nq — > for a > 0, provided that g x h 2X — g 2 h xx ^ 0. Thus, 
given Y x , the variance of Y r /d asymptotically is unbounded, 
and hence, LY/<iJ mod 2 tends to a Bernoulli 1/2 random 
variable. This proves ( fT9) . A similar analysis proves that the 
rate of the second user also improves by one bit. 
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